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Abstract 

It is widely accepted that the free energy F[H) of the two-dimensional 
Ising model in the ferromagnetic phase T< Tc has an essential branch cut sin- 
gularity at the origin H = 0. The phenomenological droplet theory predicts 
that near the cut drawn along the negative real axis H<0, the imaginary 
part of the free energy per lattice site has the form ImF[exp(±i7r) | | ] = 
±B\H\exp{—A/ I H |) for small \ H \ . We verify this prediction in ana- 
lytical perturbative transfer matrix calculations for the square lattice Ising 
model for all temperatures 0<T<Tc and arbitrary anisotropy ratio Ji/J2- 
We obtain an expression for the constant A which coincides exactly with the 
prediction of the droplet theory. For the amplitude B we obtain B = 7rM/18, 
where M is the equilibrium spontaneous magnetization. In addition we find 
discrete-lattice corrections to the above mentioned phenomenological formula 
for ImF, which oscillate in H~^. 

KEY WORDS: Ising model, droplet singularity, metastable state, false 
vacuum decay 
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I. Introduction 



It is well knowndl-i, that in the ferromagnetic phase < T< Tc the free 
energy of the two-dimensional Ising model as the function of the magnetic 
field H has a so-called droplet singularity at the origin H = 0. This sin- 
gularity prevents analytical continuation of the free energy from positive to 
negative values of H along the real if-axis. The phenomenological droplet 
(nucleation) theory^iHl^ claims, however, that the free energy can be contin- 
ued from positive to negative magnetic fields along a circle going around the 
origin in the complex if-plane (see Fig. 1). According to this theory, the free 
energy per lattice site F{H) continued in such a way gains on the negative 
real axis H < a nonzero imaginary part, which is expected to have the 
form 

ImF[exp(±i7r) \ H \ ] = ±B\H\ exp(-A/ | H |) (1) 

for small \ H \ . The sign of this imaginary part depends on the side, 
from which one approaches to the negative real axis H < 0. Expression 
([!) extrapolates to the ferromagnetic Ising model the results obtained in the 
semiclassical nucleation field theory analysis of the coarse-grained Ginzburg- 
Landau model. (0-0) In the nucleation theory, the free energy continued to 
the cut if < is interpreted as the free energy of the metastable state: 

FUH) = F{e'- \H\). 

Langer conjectured,© that lmFms{H) may be identified (up to a dynami- 
cal factor) with the metastable phase decay rate provided by the thermally 
activated nucleation of the critical droplet. 

The phenomenological droplet theory prediction for the amplitude A in 
i) is0 

A-^A (2) 
8M' ^ ^ 

where M is the spontaneous magnetization, and denotes the square of 
surface free energy of the equilibrium-shaped droplet divided by its area. 
Both and M relate to the equilibrium zero-field state, and are known 
exactly. The linear depending on prefactor in (|1|) arises in the continuum 
droplet field theory(0~B from the contribution of the surface excitations of 
the critical droplet. Voloshin claimed® that, if fluctuations are continuum 
and isotropic, the prefactor in (|l]) becomes universal. Extrapolation of the 
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Voloshin's continuum droplet field theory result to the d = 2 Ising model 
leads to the following prediction for the amplitude B : 



M 

B^By = -. (3) 

In the continuum field theory, the imaginary part of the free energy ap- 
pears in the functional integral calculations. In the alternative approach to 
the droplet singularity problem, one deals with eigenvalues of the Ising model 
transfer-matrix. Numerical transfer-matrix calculations initiated by Privman 
and Schulman©-* and continued by Giinther, Rikvold and Novotny(lil'll3) 
confirm equations (|l|) and (H) . These equations were confirmed also by Lowe 



and Wallace, 0) and by HarrisClP in numerical analysis of the small- if power 
expansion for the magnetization M{H). Recently analytic transfer-matrix 
derivation of equations (P, @ for the d = 2 Ising model has been done*-® 
in the extreme anisotropic limit. 

In this paper we generalize the transfer matrix approach developed in 



paper^c^ and verify analytically the droplet theory predictions (||), (0) for 
the square lattice Ising model for all temperatures 0<T<Tc and arbitrary 
anisotropy ratio Ji/J2- We obtain an expression for the constant A which 
coincides exactly with the prediction of the droplet theory. For the ampli- 
tude B we find B = vrM/lS, which is very close to Voloshin's result (^: 
B / By = 7r^/9 ~ 1.0966. We suppose, that this small discrepancy results 
from approximations used in our calculations. Obtained values for the am- 
plitude B are compared with those extracted numerically from the known 
coefficients of the expansion of the magnetization in powers of H by means 
of dispersion relations. 'EI ) 

We find also the discrete-lattice corrections to the phenomenological for- 
mula (|1|), which oscillate in H~^. The period of oscillations agrees well with 
that observed by Giinther et al.'^^ in numerical constrained transfer matrix 
calculations. 
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II. Transfer matrix and Hamiltonian 



The nearest neighbor Ising model on the square lattice in the magnetic 
field H is defined by the energy 



N M 



(4) 



n=l m=l 



where cr^.n = ±1, the first/second index of <Jm,n specifies the row/column of 
the lattice, Ai and J\f denote the number of rows and columns in the lattice, 
respectively. Periodic boundary conditions are implied. 

The row to row transfer matrix may be defined as T = e'^T2 Ti, where 



fi = [2sinh(2iri)]^/^exp 



^ Af \ /AT \ 

^ n=l / \ n=l / 



u 



n=l 



(5) 



Here we have used the standard notations 

Ki = /? Ji, K2 = /3 J2, h = (3H, 2Kl = -\n (tanh Ki) , 

j3 is the inverse temperature, a" (a = 1, 2, 3) are the Pauli matrices relating 
to the cite n in the row. 

The transfer matrix may be chosen in the symmetric form 



T, 



(0) 



nl/2 



T, 



(0) 



nl/2 



where T^^ is the symmetric transfer matrix of the Ising model in zero mag- 
netic field: 



T, 



(0) 



As it was shown by Schultz, Mattis and Lieb,*c3^ the latter becomes diagonal 
in fermionic variables: 



T, 



(0) 



n 



(0) 



Cexp(-7^W), 

^ ^ co{9) i;\e) m, 



(6) 



where Ti'^"-' is the zero-field Hamiltonian, 6 is the quasimomentum, C is an 
insufficient numerical factor. Fermionic operators ip^O), ip{9) satisfying the 
canonical anticommutational relations 

can be expressed in terms of the initial Pauli matrices by use of the Jordan- 
Wigner and duality transformation (see Appendix]^). The fermionic spec- 
trum uj{6) is given by 

expuj{9) = cosh 2fr * cosh 2 - cos ^sinh2ir*sinh 27^2 + (7) 
[(cosh 2K * cosh 2K 2 - cos ^ sinh 2K * sinh 2X2)^-1]^^^. 

Operator U defined by (^) also can be represented in the thermodynamic 
limit A/" ^ 00 in the fermionic variables: 

P 



where 



tr = AM^:exp^ii:, (8) 



j<n 

= i [' ^ expim f-{e) l-m + 'pH-e)] , (11) 



— 7T 



\Z2 + Z2 —2 COS 6 J 
Z\ = tanhK ^/ tanhK 2, Z2 = tanhi^^ tanhi^2; (12) 

and M is the zero-field magnetization. In the ferromagnetic phase M = 
[1 — k"^]^^^ , and k < 1, where k = (sinh 2K i sinh 2/^ 2)"^ • We have used 
in (§) the conventional notation : ... : for the normal ordering with respect 
to the fermionic operators ipiO), ip\0). Derivation of equations (^-([TT]) is 
described in Appendix Kl in the main points of which we follow Jimbo et 
alM) 

At zero magnetic field, the Hamiltonian 7i*^°^ of the Ising model is given 
by d^). Two ferromagnetic ground states | 0+) and | 0_) coexist in the ferro- 
magnetic phase k < 1. They are distinguished by the sign of the spontaneous 
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magnetization (0± | \ 0±) = ± M. The state | 0+) characterized by the 
positive magnetization +M is the ferromagnetic vacuum of '?/'(6')-operators: 

tpie) 1 0+) = for all e. 

A small magnetic field /i 7^ changes the Hamiltonian T-C^^^ to 

nih) = - In (e-^'°'/2e^ e-^'"^/^) , H(0) = H^'l (13) 
It can be expanded in powersQ of h: 

00 

n{h) = Y,'H^'\ (14) 

where Ti'--'^ ~ hK 



III. Modified perturbation theory 

Let us consider the eigenvalue problem 

n{h) \cp^{h)) = E{h)\cp^{h)), (15) 

where | 0+(O)) =| 0+). 

If /i > 0, the eigenvector | (f)_^_{h)) is the ground state of the Hamiltonian 
([T3|), and the corresponding energy E{h) is directly related with the Ising 
model free energy per lattice cite F{h,P): 

F{h,(3)=F{0,f3) + ^. (16) 

The energy can be expanded in the formal power series 

00 

i=i 

which coefficients Cj can be, in principal, determined from standard Rayleigh- 
Schrodinger perturbation theory. 

^In the extreme anisotropic limit^l^^ the Hamiltonian expansion containes only 
two terms: n{h) = n^°^ + H^^', where H^^) = -U. 
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However, if the magnetic field is small and negative h < 0, the state 
I 4>+{h)) (with positive magnetization almost equal to M) must be identified 
with the metastable (false) vacuum. It decays due to the quantum tunneling, 
and the decay rateQ F is proportional to the imaginary part of the energy 
E{h) continued to negative magnetic fields: ^'^3^ 

r = -2 Im E{h). (17) 

It turns out, however, that F can not be determined from the straightforward 
perturbation theory with the zero-order Hamiltonian 'H^'^\ This is due to 
the fact that the term Ti.^^^ in the expansion ([T^) contains the long-range 
interaction {—Uq) between fermions, which is given by 



-Uo = -U Uie)^H I M J] : exp -2 5^ b]b, : . (18) 

ngZ \ j< n ) 

Here 

^^\e) exp i-zjo) , 6, = y -m exp m , 

are the operators which create/annihilate a fermion at the cite j. Interaction 
( [TE| ) increases linearly with the distance between fermions,®^ and therefore 
changes the structure of the Hamiltonian spectrum. So, to describe decay of 
metastable vacuum | 0_,_(/i)), one should include the long-range interaction 
([TSp into the zero-order Hamiltonian. 

Accordingly, we subdivide the Hamiltonian (|1^ into the zero-order TYq 
and interaction V parts, as follows: 

n{h) = no + v, (19) 

where 

Ho = -f/o-AT \h\ M, (20) 

CO 

V = n'^^^ + Uo+Af \h\ M + (21) 

i=2 

The numerical constant JV \h\ M in ( PD| ) is chosen to provide TYq | 0+) = 0. 

^Strictly speaking, the term "decay rate" here relates to the quantum-mechanical model 
with Hamiltonian (jlj), but not to the initial two-dimensional Ising model (jj). The nucle- 
ation rate in the latter model contains also the so-called kinetic prefactor, which depends 



on the detailed non-equilibrium dynamics.'!^' 
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A. Zero order spectrum 



Consider the zero-order eigenvalue problem 

Ho\ct>i) = Ei \ct>i) (22) 

First note, that eigenstates | (pi) can be classified by the fermion number, 
since the modified zero-order Hamiltonian ( PP] ) conserves the number of 
fermions. As in paper, ^IHI) we shall consider only two-fermions (i.e. one- 
domain) states in (0). Physically, this means that we neglect interaction 
between nucleating droplets of the stable phase. 

In the coordinate representation equation (^) takes the form 

^ Knn' <Pi{n') - M\nh \ (pi{n) = -^(pi{n), 

n'GZ 

where 

<Pi{n) = (0+ I ^0 &n I <Pi{-n) = -(piin). 

If the energy Ei is small enough Ei <^ ^^(0), the wavefunction 0;(n) is 
mainly concentrated far from the origin in the classically available region 
I n \> a;(0)/(| h \ M). Therefore, we can apply the 'strong coupling 
approximation' to represent the wavefunction in the form 

0,(n) ~ (/?;(n) -^;(-n), (23) 

where the function ifiin) solves the equation 

^ Knn' \ h\M n (pi{n) = Y<^/(n). 

n'eZ 



After the Fourier transform, we obtain 



= I —(pi{e) exp (ine) , 
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where 



cp,{e) = c exp^-^^j-yj^[fie)-Eie]\, (24) 

C = {2\h\M Afy^/^ 
f{e) = 2 f dx^ix)- 



The 27r-periodicity condition for the function ^Pi{9) determines the energy 
levels El. 

El = -2\h \ ML (25) 

TT 

The normalization constant C in (^) is chosen to yield 

where AE = 2 \ h \ M is the interlevel distance. 

B. Decay rate 

The first-order correction to the false vacuum energy is trivial: E'^^^ = 
(0+ I V I 0+) = Af \h\ M, the second-order correction is given by 

E(^) = -Ai^5^^^L^^^ (26) 

To determine the decay rate of the false vacuum, the following trick is used. 
We shift the excitation energy levels Ei in (^) downwards into the complex 
i?-plane: Ei ^ Ei — i'-f, where the width 7 describes phenomenologically the 
decay rate of one-domain states | 0^). Decay of these states should be caused 
by the interaction term (pT]) in the same manner as the false vacuum decay . 
As the result, the metastable vacuum energy gains the imaginary part 

Im E ~ -TT g{h) \{<Pi\V\ 0+) |%=o, (27) 

where 



g{h) = Imcot 



/(tt) - iivj 
2 \h \ M 



(28) 
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The metastable vacuum relaxation rate T is determined then in the usual way 
([T7|). It is evident from (pTj) , (|1^ that F oscillates in with the period 
Ah^^ given by 

Ah-^ = 27rM//(7r). (29) 

These oscillations become considerable in the case of the resonant tunneling 
7 < AE. In the opposite limit 7 ^ AE oscillations in vanish and 
relations (p?]), ( p!?] ) transform to the Fermi's golden rule:*!^'^^ 



27T 



E,=0 ■ 



(30) 



Let us now calculate the matrix element in (pT]): 



V\0+) 



de 
2^ 



<pm (o+i^(-^)^(^^) v|o+) 



-(^r(^)(o+ I ^(-e)^(^)i^|o+). 

— TT 

Expanding operator V in the /i-power series 



(31) 



(0+ 1 ij{-9) m y 1 0+) = 5^(0+ 1 ^(-^) m n^^^ 1 0+) 

and keeping in it only the leading (j = 1) term one obtains from (P), (13), 



(Oh 
(Oh 



i){~6) ij{9) 1^ I 0+) ~ (0+ I ^{-9) ij{9) I 0^) 



ij{-9)mu\0+) 



u;{9) 



sinhu;(^^) 



-2W\h \ M 



dlne{9) uj{9) 



d9 sinha;(^) 



Thus, the matrix element (^) can be approximately represented as 



h I V I 0+) 



iAf\h\M 



d9 u{9) d\nei9) 



TT sinhu;(6') 



d9 



(32) 



Substitution of (|32D , (|27|) , (^4|) into (|T6D yields the final expression for the 
imaginary part of the free energy F^s in the limit h —0: 

2 



±111 J. — 



TT 



H I M g{h) 



d9 oj{9) d\ne{9) 
TT sinhco'(6') d9 



exp 



2 U I M 



(33) 
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This expression generahzes formula (18) of reference^H* to arbitrary anisotropy 
Ji/ J2 and all temperatures < T < Tc. 

The last (exponent) factor of the integrand in (|33|) oscillates with high 
frequency in the considered case of small \h\. Therefore, in the limit \h\ — >■ 0, 
the integral in the right-hand side of (|33D is determined by the saddle point 
of fiP): 9 = 9i = —ilnzi, u!{9i) = 0, and asymptotically equals to 

B\H\ g{h) exp 



Im F^s ~ 



H 



where 



A 
B 



M 13 
ttM 



(34) 

(35) 
(36) 



IV. Discussion 



First, let us establish equivalence of expressions (g) and (pSQ for the am- 
plitude A, which are given by the phenomenological droplet theory and by 
our transfer-matrix calculations. 

The droplet equilibrium shape in the d = 2 Ising model is determined by 
the equation 

ai cosh(/?Axi) + 02 cosh(/?Ax2) = 1, (37) 

obtained by Zia and Avron/0) Here Xi,X2 denote Descartes coordinates of a 
point on the droplet boundary, the scale parameter A determines the droplet 
size, and 



ai 



tanh(2is:2) 



a2 



tanh(2iri 



cosh(2iri) ' cosh(2is:2)' 
It is remarkable, that equation (^) can be rewritten in terms of the Ising 
model excitation spectrum (0) first obtained by Onsager:*^© 



X2 = ±— uj{i(3Xxi). 



(38) 



Integrating in Xi this equation we find the area of the equilibrium-shaped 
droplet 5(A) = W/X^, where 

2 

J' 



W 
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It follows from Wulff's theorem'^EZl* that the surface energy S(A) also can be 
expressed in W: S(A) = 2W/\. Therefore, = AW, and 



8 M M 13 

in exact agreement with (pSf). 

Our expression vrM/lS for the amplitude B is the same as that obtained 
previously in the extreme anisotropic limit. As it was mentioned in the 
introduction, this expression is very close to the Voloshin's result (^. The 
latter is expected to be exact in the critical region, where fluctuations are 
isotropic and universal. It is likely, that the small discrepancy between (|36D 
and @ is caused by approximations used in our modified perturbation theory. 
We hope to clarify this question in future. 

Let us compare obtained expressions for the amplitude B with the numer- 



ical results by Baker and KimM which they calculated for the symmetric 
case Ji = J2 = J of the Ising model . These authors considered the power 
series for the magnetization M{h) at a fixed temperature below T^: 



M{h) = M(0) -2^(-2/i) 



n 



n=l 



CI 



and calculated numerically 12 coefficients a„ in this expansion at m = 0.1 u 
and 24 coefficients at u = 0.9 Uc- Here u = exp(— 4/5J); Uc = S — -\/8 corre- 
sponds to the critical temperature. On the other hand, Lowe and Wallace© 
demonstrated by use of the dispersion relation, that equation (|l|) leads to the 
following asymptotic formula for the coefficients a„: 

an - ^(2A/3)-"^^^±^. (39) 
n^oo 2tt n 

So, the ratio 

R - ^(^ + 1)! un) 

" " 27rn a„ (2A/5)- 

should approach to unity at large n, if we put in it the correct values of 
amplitudes A and B . We plot in Fig. 2 this ratio, where coefficients a„ were 
taken from paper*^© by Baker and Kim. The left pair of curves corresponds 
to the low temperature case u = 0.1 Uc, the right pair of curves corresponds 
to the higher temperature u = 0.9 Uc- The amplitude A in (HO) is taken from 
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Solid and dashed curves differ by choice of the amphtude B in 
In sohd curves, it is chosen as B = ttM/18 according to our resuh (|3^); in 



dashed curves B = M/ (27r) according to Voloshin's result. ^lLhP All four curves 
in Fig. 2 seems to stabilize at large n to the values, which are rather close to 
unity. This indicates a remarkable good agreement of numerical results^cil^ 
with expressions or (^). Though, agreement with Voloshin's value seems 
somewhat better, saturation in curves is not achieved, and further numerical 



calculation are desirable to distinguish between and (|35|). 

Expression (|3^) differs from (|l]) by the oscillating factor g{h). We inter- 
pret this factor as the correction coursed by the discrete-lattice effects. These 
oscillation being negligible in the critical region, may be significant at low 
temperatures, especially in the presence of strong anisotropy.^^^l) The period 
of oscillations in is given by (pOf). It is plotted in Fig. 3 in the symmetric 
case Ki = K2 = K. Such oscillations with period Ah~^ 1/2 were observed 
in numerical constrained transfer matrix calculations by Giinther et al.^-^ 
aX K = 1. This period agrees well with our value A/i^^ = 0.494891, which 
follows from (120). 
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Appendix A. Fermionization 

In this Appendix we present fermionic representations of spin operators, 
which are used in Section |T|. Consideration is restricted to the ferromagnetic 
phase T < Tc in the thermodynamic hmit A/" — oo. In this limit the Jordan- 
Wigner transformation can be written as^"' 



P. 



(Al) 



Here P„ and are the fermionic operators satisfying the following anticom- 
mutational relations: 



{Pn , -Pn'} — 25^ 



{Qn , Qn'} — — 25„ 



{P„,Q„,} = 0. 



Let us define the another set of fermionic operators p„, g„, which are related 
with Pn , Qn by the duality transformation:^® 



Pn 



^ Qni Qn 



-iP 



n+l- 



(A2) 



Operators p„, qn obey the same anticommutational relations as Pn , Qn, 
span an orthogonal space of free fermion field, and generate the Clifford 
algebra*!^ . 

Fermionic creation and annihilation operators ip{9),ip\0) introduced in 
Section |T| are related with p„, g„ by 



where 



p{0) 



(A3) 



'ind 



Pn 



9W = E^"'"'^" 



— tanh K2 



-1, 



zi){e 



id 



Z2){e 



ie 



1/2 



(A4) 



parameters Zi, Z2 are defined by (|T^) . 

Relations (|A1|) , (|A2|) express operators p„, qn in terms of Pauli matrices. 
The inverse transformation reads as 



0". 



Pn Qn-l, 

Pn {Pn-1 qn-2) {Pn-2 Qn-s) ■ ■ ■ 
iqn-1 {Pn-1 qn-2) {Pn-2 qn-s) 



(A5) 
(A6) 
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Brackets in equations (^), ( |A6| ) are shown to indicate, that cr^ and are 
the products of odd number of fermionic operators, i.e. and are the 
odd elements of the Chfford group. 



Following Jimbo et al.^^ , let us introduce operator a^, which represent 
cr^ under the boundary condition cr^ — ^ 1 for n — — cx): 

al = Pn-l) {qn-2Pn-2) ■■■■ 

This operator is an even element of the Clifford group. Due to the obvious 
identity cr^cr^, = a^cr^, , operators a\ and a\ produce the same correlation 
functions, which makes reasonable to identify them in the thermodynamic 
limit. Accordingly, we shall replace cr^ by a\ in the operator U : 

U^hJ2^l (A7) 

n=l 

Operators cr^ are characterized up to ±1 factor by the following commu- 
tation relations: 

al qn' = x{n - n) qn' , p„/ = x(n - n') (A8a) 

where 

, . f 1 for n > 
^ ' [ - 1 for n < 

Thus, induces a linear orthogonal transformation of the linear space of 



free fermions. As it was shown by Jimbo et al.^^ in Appendix 1, such an 
operator can be expressed as the normally ordered exponent 

a^ = (a3):exp(p„/2):. (A9) 

Here operator p„ is quadratic in free fermionic variables p„ and g„, (cx^) 
is the vacuum expectation value of d"^, i.e. the spontaneous magnetization 
(a^) = M = [1 — k"^^^^ . The explicit expression for p„ reads as 



^ = (AlO) 



3<n 



" de 



= r ^ e^^' U49) qie), 
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where the functions 



provide the Wiener-Hopf factorization of exp [2ia 

exp[2ia(^)] = U+{e) U.{9) 



Functions U+{9) and U_{d) are analytical m. z = e*^ outside and inside the 
unit circle, respectively. Rewriting ip'^^'^ and ■* in terms of creation and 



annihilation operators V'^(^); V'(^) by use of ( [A^ ) one obtains from 
(|A9| ) , ( |A10| ) the desired fermionic representation (|^) of the U operator. 

In deriving ( |A10D we have chosen the free fermion basis p(^), q{0) and 
applied the theorem, presented by Jimbo et al. in pages 137, 138 of their 
article^c^l) . In our case, the kernel functions for matrices P, E introduced 
in this theorem read as (compare with equations (3.15), (3.16) in the same 
article): 



gi(n-i)(e-e') 

X _ g-t(6»-e'-io) ' 

exp 
exp [2ia{ 







2n5{e - e') 



To illustrate convenience of representations (H), (|A9|) , (|A10|) , we shall 
apply them to derive a compact Fredholm determinant formula for the zero- 
field correlation function (o"o,o c"m,n) in the ferromagnetic phase. Let us first 
write this correlation function by use of ( |A9| ) , (|^) in the form 



(0+ 



al exp( 



T^W) al I 0+) 



(Oh 

: exp (po/2) : exp (-T^^m) : exp (p„/2) :| 0+), (All) 



where 



Po 
2 

Pn 

2 



^ {2ixf ^ 



{ei.e^) ^(^1)^(^2) 



^ dOide^ 
{2Txf 

1 



MM+e,) ^-in(B,+e,) jj^o^^Q^) ^t(^^) ^t(^2 



2 sin [(^1 + ^2) /2] 
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We have dropped all creation operators ijj^ in and all annihilation oper- 
ators ip in since the normally ordered exponents of p and Pq act in ( All 



on the vacuum states. In the well-known holomorphic representation^&lP of 
fermionic operators, the matrix element ( [All| ) takes the form of the Gaus- 
sian continual integral over Grassnann variables, which integration yields 
immediately: 

{crofl Om,n) = det (1 - Drr^n) = (A12) 

^'-p(-E^) (A13) 

Here Dmn denotes a linear integral operator acting on the function f{9) as 
follows: 



Equation ( |A12| ) is a compact forms of the well-known exact representation 



of the two-point correlation function in the Ising model obtained by Wu et 
al.^^^ (see equations (2.9)-(2.13) in the referred article). The latter repre- 
sentation can be reduced to (|A13|) by explicit integration in the right-hand 
side of equation (2.12) in (pi, 03, (j)^ .... 



17 



References 

^N.S. Isakov, Commun. Math. Phys. 95:427 (1984). 

^D.B. Abraham and P.J. Upton, Phys. Rev. Lett. 69:225 (1992). 

^D.B. Abraham and P.J. Upton, Phys. Rev. Lett. 70:1567 (1993). 

^A.F. Andreev, Sov. Phys. JETP 18:1415 (1964). 

^M.E. Fisher, Physics 3:255 (1967). 

^H. Kastrup, Phys. Rev. Lett. 81:2280 (1998). 

^J.S. Langer, Ann. Phys. (N.Y.) 41:108 (1967). 

«N. J. Giinther, D.A. Nicole and D.J. Wallace, J. Phys. A 13:1755 (1980). 

9R.K.P. Zia and D.J. Wallace, Phys. Rev. B 31:1624 (1985). 

^oM.B. Voloshin, I.Yu. Kobzarev, and L.B. Okun', Yad. Fiz. 20:1229 (1974) 
[Sov. J. Nucl. Phys. 20:644 (1975)]. 

^^S. Coleman, Phys. Rev. D 15:2929 (1977), Erratum: Phys. Rev. D 16:1248 
(1977). 

^^C.G. Callan and S. Coleman, Phys. Rev. D 16:1762 (1977). 

^^P.A. Rikvold and B.M. Gorman, in Annual Review of Computational 
Physics /, edited by D. Stauffer (World Scientific, Singapore, 1994). 

^^M.J. Lowe and D.J. Wallace, J. Phys. A 13:L381 (1980). 

i^C.K. Harris, J. Phys. A 17:L143 (1984). 

i^C.C.A Giinther, P.A. Rikvold, and M.A. Novotny, Phys. Rev. Lett. 71:3898 

(1993) . 

i^C.C.A Giinther, P.A. Rikvold, and M.A. Novotny, Phystca A 212:194 

(1994) . 

^^M.B. Voloshin, Yad. Fiz. 42:1017 (1985) [Sov. J. Nucl. Phys. 42:644 
(1985)]. 



18 



i^V. Privman and L.S. Schulman, J. Stat. Phys. 29:205 (1982). 

20s. B. Rutkevich, Phys. Rev. B 60:14525 (1999). 

21G.A. Baker Jr. and D. Kim, J. Phys. A 13:L103 (1980). 

Schultz, D.C. Mattis and E.H. Lieb, Rev. Mod. Phys. 36:856 (1964). 

23m. Jimbo, T. Miwa, Y. Mori and M. Sato, Physica D 1:80 (1980). 

^'^L.D. Landau and E.M. Lifshitz, Course of Theoretical Physics, (Pergamon 
Press, Oxford, 1965), Vol. 3. 

2^M. Reed and B. Symon, Methods of Mathematical Physics, (Academic, New 
York, 1978), Vol. 4. 

2^J.M. Ziman, Principles of the Theory of Solids (University Press, Cam- 
bridge, 1972). 

2^R.K.P. Zia and J.E. Avron, Phys. Rev. B 25:2042 (1982). 

Onsager, Phys. Rev. 65:117 (1944). 
29R. Kotecky and E. Olivieri, J. Stat. Phys. 70:1121 (1993). 
=^°J.B. Kogut, Rev. Mod. Phys. 51:659(1979). 

3^A.A. Slavnov and L.D. Faddeev, Introduction to Quantum theory of Guage 
Fields (Nauka, Moskow, 1988). 

32T.T. Wu, B.M. McCoy, C.A. Tracy and E. Barouch, Phys. Rev. B 13:316 
(1976). 



19 



Figure Captions 



Figure 1: Free energy continuation paths from the positive real 
axis H > to the cut going along the negative real axis H < 0. 

Figure 2: Plot of Rn given by (^) versus n. Coefficients a„ in 
(^) are taken from reference amplitude A is taken from (p5|); 
amplitude B is taken either from (^) (solid curves), or from (|]) 
(dashed curves). Two left curves correspond to u = O.lUc, two 
right curves correspond to u = 0.9uc- 

Figure 3: Oscillation period Ah^^ versus K in the symmetric 
case. 
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